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ABSTRACT 

Second-order  steady  wave  forces  on  a  restrained,  submerged, 
circular  cylinder  are  studied  experimentally  as  a  function  of 
cylinder  depth  and  wave  frequency.  The  experimental  results  are 
compared  to  forces  calculated  from  a  linearized  potential  theory. 
It  is  found  that  in  all  cases  the  experimentally  observed  forces 
are  higher  -  in  some  cases  by  an  order  of  magnitude  -  and  are  not 
attenuated  with  depth  as  rapidly  as  predicted  by  the  theory. 
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1.  INTRODUCTION 

Since  very  early  in  written  history  warriors  have 
attempted  to  escape  the  sea's  surface  and  utilize  submarine 
warfare.  Herodotus  (460  B.C.),  Aristotle  (332  B.C.),  and 
Pliny,  the  elder  (77  A.D.)  mention  determined  attempts  to 
build  submersibles .  Alexander  the  Great  (35&  to  323  B.C.) 

is  the  first  person  known  to  have  descended  into  the  sea  in 

1 

a  vessel  of  any  kind. 

Since  that  time  submarines  have  been  invented  which 
were  propelled  by  oars ,  sails ,  treadles ,  hand-operated  screws , 
clockwork,  springs,  steam  stored  in  tubes,  chemical  engines, 
compressed  air,  stored  gases,  electric  motors,  and  nuclear 
power.  Over  three  hundred  years  ago  Mother  Shipton,  a  famous 
English  prophetess,  predicted  the  coming  of  the  submarine  when 
writing,  "underwater  man  shall  walk,  shall  ride,  shall  sleep, 
shall  talk."  Her  prophesy  has  come  true,  but  in  carrying 
out  her  prophesy  man  has  discovered  that  even  though  he  may 
be  an  underwater  man  in  many  senses,  he  has  not  completely 
escaped  the  forces  of  surface  waves. 

Man  was  early  to  recognize  the  forces  of  the  sea  to  be 
powerful,  at  times  destructive,  and  inevitable.  He  set  out  to 
study  them,  attempt  to  predict  them,  and  to  design  vessels 
which  would  withstand  and  take  advantage  of  them.  Studies  of 
forces  on  vessels  penetrating  a  free  surface  are  too  numerous 
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to  mention  here. 

The  United  States  Navy  commissioned  the  USS  Holland,  its 
first  submarine  on  April  11,  1900,  The  military  importance  of 
submarines  has  increased  at  a  staggering  rate  since  that  time. 
More  recently,  commercial  and  scientific  interests  have  awakened 
to  the  vast  resources  under  the  sea,  hence  the  use  of  non- 
military  submersibles  has  begun  to  increase  at  a  rate  that  may 
exceed  that  of  its  military  cousin.  As  more  and  more  use  is 
made  of  undersea  vessels,  they  become  larger  —  the  largest 
known  maneuver  able  undersea  vessel  to  date  is  about  ten  meters 
in  diameter,  greater  than  one  hundred  meters  in  length,  and 
displaces  in  excess  of  six  thousand  tons  —  and  yet  they  must 
be  precisely  controlled. 

In  order  to  be  able  to  provide  this  precise  control,  the 
forces  acting  on  these  vessels  must  be  known,  One  source  of 
such  forces  is  surface  waves.  A  number  of  studies  have  been 
made  of  the  wave  induced  forces  on  underwater  bodies ,  but  most 
of  these  "studies  have  been  concerned  with  an  underwater  body 

moving  with  some  translational  velocity  with  respect  to  the 

2 

fluid.  Examples  of  these  studies  are  some  of  Havelock's 

?  U> 

papers,  Kim's  J   paper,  and  Salvesen's  paper. 

The  present  work  will  concern  itself  with  a  body  that  is 
not  moving,  but  which  is  rigidly  restrained.  Kulin  studied  the 
wave  forces  on  submerged  cylinders  and  plates  by  considering 
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the  effect  of  a  single  intumescence  which  subjects,  the  fluid  to 
acceleration  followed  by  deceleration.  He  studied  (l)  the 
extent  local  acceleration  causes  the  drag  coefficient  to  depart 
from  steady  state  values,  (2)  the  extent  that  friction  affects 
potential  flow  values  of  inertial  coefficient,  and  (3)  the  extent 
of  force  dependence  upon  preceding  flow  history,  as  well  as  on 
instantaneous  conditions.  He  found  that  the  forces  could  not  be 
represented  by  a  single  universal  total  resistance  coefficient, 
and  that  the  ratio  of  the  distance  of  fluid  motion  to  object 
size  was  a  significant  parameter  related  to  vortex  formation 
during  fluid  deceleration. 

In  19&2,  Ogilvie  applied  and  extended  the  form  of 

7 
solution  that  had  been  developed  by  Ursell  in  1950  to  calculate 

the  first-order  oscillatory  force  and  the  second-order  steady 
force  in  the  following  situations:  (a)  a  submerged  cylinder  is 
restrained  from  moving  under  the  effect  of  incident  sinusoidal 
waves;  (b)  a  submerged  cylinder  is  forced  to  oscillate  sinusoidally, 
in  otherwise  calm  water;  (c)  a  cylinder,  which  is  neutrally 
buoyant,  is  allowed  to  respond  to  first-order  oscillatory  forces 0 
He  proves  that  subject  to  his  assumptions  a  knowledge  of  the  first- 
order  potential  supplies  information  sufficient  to  solve  these 
problems.  This  paper  deals  with  an  attempt  to  experimentally 
verify  his  results  for  case  (a),  the  restained  cylinder,, 

Assume  that  a  circular  cylinder  is  located  under  a  free  surface, 
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(see  figure  1.)  such  that  the  surface  of  the  cylinder  is 
specified  by 

S(x,y)  =  x2  +  (y  +  h?    -  a2  =0„  (1) 

The  undisturbed  free  surface  is  taken  as  the  x-axis  and 
the  instantaneous  free  surface  10.11  be  specified  by 

y  -  Y(x,t)  =  O0  (2) 

The  required  velocity  potential,  JD(x,y,t)  ,  satisfies 

\*  +  V  =  °  for     y  <  Y(x»t),  (3) 

and       S(x,y)  >  0, 
!D/YX   -  1DY  +  Y(    =  0  on       y  =  Y(x,t)  ,  (4) 

gY  +  Et  +  i(ffij  +  ffl})  =  0       on       y  =  Y(x,t) ,  (5) 


and 


EXSX   +  EySY  +  St   =  0       on       S(x,y)   =0,  (6) 


plus  appropriate  conditions  at  infinity.  The  time  average,  Zu,(t)  , 
of  the  second-order  steady  force  can  be  shown  to  be 


zPTty*=  stej*-  iY^cty* 


d©  (7) 


or  if  we  let  55(x,y,t)  be  the  real  part  of  a  function  of 


-8- 


a  complex  variable,  f(z,t),  where  z  =  x  +  iy,  then 

t      r  j 

dO  (8) 


¥%£)*=  -%1a.fi  ei0  [f '  ( z ,  tjf^.tjj^c 


.7T 


=  -i27Tps(KO)  e     Q(ka,kh).  (8«) 


k  is  the  wave  number  defined  as 

k  =^/g,  (9) 

and  HO  is  the  amplitude  of  the  incident  waves. 

Define  steady  vertical  force  coefficient ,p  ,   as 
.t 


B  =  5551      -  do) 

2irfg(mf 


From  (8*)  and  (10),  it  is  seen  that 
-2kh 


j3  =  e~^   Q(ka,kh).  (11) 

It  is  this  quantity  obtained  by  Ogilvie  that  will  be 
compared  with  experimental  results  in  Section  IV. 


/ 
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Figure  1.  Geometry  of  the  probl 
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II.  TEST  APPARATUS  AND  PROCEDURE 

II. A.  TEST  APPARATUS  DESCRIPTION 

This  experiment  was  conducted  in  the  wave  tank  in  the  M.I.T. 
Department  of  Naval  Architecture  and  Farine  Engineering  Hydro- 
dynamics Laboratory.  The  working  section  of  the  tank  is  ten 
feet  long,  one  foot  wide,  and  six  inches  deep.  The  wave  generator 
is  a  vertically  oriented  pusher,  hinged  at  the  bottom,  and  driven 
by  a  push-rod  connected  off  center  on  a  flywheel,  which  is  in  turn 
driven  by  a  variable  speed  electric  motor.   The  observation  section 
was  located  about  equidistant  from  the  two  ends.  An  energy 
absorbing  beach  made  of  rubberized  horse  hair  was  located  at 
the  end  opposite  the  wave  generator.  This  beach  was  not  completely 
successful  in  preventing  reflected  waves.  This  will  be  considered 
in  Section  II. 3. 

Wave  height  was  measured  by  a  capac: tance  probe  mounted  on 
a  carriage  which  was  driven  by  a  variable  speed  electric  motor 
and  ran  on  tracks  along  the  top  of  the  tank.   The  output  of  the 
probe  was  fed  through  a  buffer  amplifier  to  a  Hewlett  Packard 
Model  7702  B  recorder.  The  carriage  was  moved  at  a  constant 
speed  before  each  force  observation  to  obtain  data  for  cal- 
culating the  reflected  to  direct  wave  ratio  and  then  placed 
at  a  predesignated  point  as  close  as  possible  to  the  cylinder 
to  measure  wave  height  during  the  period  of  force  measurement. 
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The  cylinder  itself,  a  one  inch  diameter  lucite  rod  11  15/l6 
inches  long,  was  suspended  below  the  surface  by  a  3/32  inch 
diameter  rod  tapped  into  the  cylinder  at  points  1/8  inch  from 
each  cylinder  end.  Any  effect  on  vertical  forces  by  these  thin 
supporting  rods  was  assumed  to  be  negligible.  The  l/32  inch 
clearance  between  each  end  of  the  rod  and  the  tank  wall  is 
considered  to  be  smaller  than  any  wall  boundary  layer  thickness, 
and  it  is  felt  that  there  was  no  appreciable  departure  from  the 
two  dimensional  condition.  The  tops  of  the  supporting  rods 
were  threaded  and  attached  to  a  yoke  which  was  restrained  so 
that  its  only  degree  of  freedom  was  in  the  vertical  direction. 

Vertical  force  on  the  yoke  was  measured  by  two  Schaevitz 
Engineering  Kodel  No.  FTA  3  Force  Transducers,  one  on  each  end. 
The  electrical  outputs  of  these  transducers  were  connected  and 
processed  as  shown  schematically  in  figure  2.  As  seen  from  the 
figure  the  output  voltage  was  directly  proportional  to  the  time 
integral  of  the  steady  force  component,  a. 
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II. B.  PROCEDURE 

As  mentioned  in  Section  II. A.,  the  energy  absorbing  beach 
was  not  completely  successful  in  preventing  reflected  waves  and 
it  was  considered  necessary  to  take  these  reflected  waves  into 
account.  A  method  of  doing  this  by  the  use  of  a  standing  wave 
ratio,  p,  is  developed  in  Appendix  A.   That  result  is  given  in 
(A7)  and  is- repeated  here, 

j9=-j£° r  (12) 

The  values  of  HX  and  Wi   used  in  (A8)  were  obtained  by 
allowing  the  capacitance  probe  carriage  to  move  at  constant 
speed  along  the  tracks  as  mentioned  in  Section  II. A.  and  observ- 
ing the  maximum  and  minimum  values  of  wave  height  recorded. 

The  wave  height  recorder  was  calibrated  by  changing  the 
still  water  level  by  a  known  amount  and  observing  the  change 
in  recorder  level.  Because  of  the  sensitivity  of  a  capacitance 
probe,  the  calibration  was  repeated  frequently  and  the  appropriate 
calibration  factor  was  used  in  each  calculation  of  A 

Force  measurement  calibration  was  accomplished  by  first 
opening  the  shorting  switch  shown  bypassing  the  integrating 
amplifier  in  figure  2.  for  a  known  period  of  time,  TD,  and 
observing  the  integrated  drift  voltage,  ED.   The  switch  was  then 
closed  and  a  known  weight  Y  was  ^laced  on  the  cylinder  snTvnoT,tincr 
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yoke.  The  switch  was  reopened  for  a  known  period  of  time,  TS , 
and  the  integrated  signal  voltage,  ES,  was  observed.   The  cali- 
bration factor,  AE,  was  then  determined  by  the  following  formula: 

AE  =  Y/E  (13) 

where 

E  =  (ES/TS)  -  (ED/TD).  (1*0 

This  factor  was  found  to  remain  constant  at  O0O36  pounds  per  volt. 

The  unknown  steady  wave  force,  Y^Kt) ,  was  found  by  the 
following  equation,  using  values  of  ES,  TS,  ED,  TD,  HM,  and  HX 
observed  for  various  cylinder  depths,  h,  and  incident  wave 
frequencies,  F, 

T*TO"*=  (E)(AE)  (15) 

y#0was   calculated  using   (10)    and  (15)    and  ft  was   calculated  by  (12) 

Calculations  were  performed  at  the  M.I.T.  Computer  Center 
using  the  program  in  Appendix  B, 
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III.  RESULTS 

Experimental  data  is  tabulated  in  Appendix  C„  Values 
of  j5  calculated  from  this  data  using  (12)  are  tabulated  on 
the  following  page. 

These  values  of  /S  are  presented  graphically  in  figures 
3.  through  9.  along  with  curves  for  ka  =  0.2  and  0.5  from 
Ogilvie's  theory. 

In  drawing  the  experimental  0   curves ,  curves  for  the  next 
higher  and  lower  values  of  ka  were  considered  as  well  as  $ 
values  for  the  curve  being  drawn.  As  a  result  the  experimental 
curves  for  some  ka  values  do  not  pass  through  all  data  points. 

The  individual^  curves  from  figures  J»   through  9.  are 
reproduced  in  figure  10,  for  ease  of  comparison. 
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Figure  3.  Steady  vertical  force  coefficient,  ka  =  0.25 
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Figure  5.  Steady  vertical  force  coefficient,  ka  =  0.35 
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Figure  6.  Steady  vertical  force  coefficient,  ka  =  0.4 
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Figure  7.  Steady  vertical  force  coefficient,  ka  =  0.4-5 
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Figure  8.  Steady  vertical  force  coefficient,  ka  =  0.5 
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Figure  9.  Steady  vertical  force  coefficient,  ka  =  0.55 
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Figure  10.  Steady  vertical  force  coefficients 
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,  IV.  DISCUSSION  OF  RESULTS 

Reference  to  figures  3»  through  9»  shows  that  the 
experimental  curves  of  steady  vertical  force  coefficient  are 
in  all  cases  higher  —  in  some  cases  more  than  an  order  of 
magnitude  —  than  the  theoretical  curves ,  and  in  all  cases 
except  for  ka  =  0.25  and  0.3  the  slope  is  less  steep  than  the 
theoretical  curves. 

It  is  also  noted  that  the  experimental  curves  all  have 
inflection  points  between  the  third  and  fourth  data  points  — 
at  an  h/a  ratio  of  about  H-,     From  figure  10.  it  is  easily- 
seen  that  this  inflection  is  much  more  pronounced  in  the  curves 
for  smaller  values  of  ka  —  longer  wave  lengths.   Since  for  the 
longest  wavelength  investigated  the  tank  depth  was  over  a  half 
wave  length  it  was  expected  that  it  would  have  an  infinite 
effective  depth;  however,  a  finite  depth  effect  seems  to  be  a 
plausible  explanation  for  these  curve  inflections. 
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V.      C01ICLUSIQKS  AHD  RECOU-EMDATIOHS 

From  Sections  III  and  IV  it  may  be  concluded  that  Ogilvie's 
potential  theory  does  predict  the  general  behavior  of  the  steady 
vertical  force  coefficient,  but  that  in  all  cases  it  predicts 
a  force  smaller  than  that  measured  experimentally  and  that  the 
forces  will  be  attenuated  more  rapidly  with  depth  than  that 
experimentally  observed. 

Because  of  the  curve  inflections  noted  in  Section  IV  and 
since  smaller  values  of  ka  are  more  likely  to  be  encountered  at 
sea,  i.e.,  a  one  hundred  meter  wave  length  and  a  ten  meter 
diameter  cylinder  give  a  ka  of  0.31^  —  most  waves  of  interest 
are  longer  and  most  diameters  are  smaller,  hence  most  ka's  of 
interest  are  smaller,  it  is  felt  that  an  extension  of  this 
experiment  .to  smaller  values  of  ka  in  facility  with  an  infinite 
effective  depth  and  no  reflected  waves  would  be  rewarding. 
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APPENDIX  A 


Theory  Modification  to  Account  for  Reflected  Waves 

o 

Using  the  'circle  theorem'  of  Milne-Thompson,  the  complex 

potential  function,  f(z,t),  in  (8)  can  be  represented  by 

-fert  +  kz)       -kh  iTcrt  +       ka* 
3  +  e        e  1  (z  +  ih). 


f,  (z,t)   =  A 

where 

A  =  (HO)g/cr. 

Using  (Al) ,    (8)    and  (10),     is  found  to  be 

-2kh 
fi=  e  (ka)   I,  (2ka). 

The  complex  potential  function  for  a  reflected  wave 

traveling  in  the  opposite  direction  with  an  unknown  phase 

angle ,  H',  and  a  magnitude  of  p  times  the  direct  wave  can  be 
given  by 


(Al) 


(A2) 


(A3) 


f-  (  z ,  t)  =  pAe 


if 


-i(<rt  -  kz)   -kh  i 
+  e   e 


<rt   -   ka* 


(z  -  ih) 


(A^) 


Now  the  complex  potential  function  for  a  direct  and  reflected 
wave  existing  simultaneously  is  given  by 

f(z,t)  =  f,  (z,t)  +  fa(z,t).  (A5) 

Using  (A5),    (8)   and  (10),  00  is  found  to  be 

-2kh 


00  =  (1  +  p2)    e         (ka)   I,  (2ka). 


(A6) 
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Consequently 

0   =P0/(1  +  p2)  (A?) 

where  B  would  be  the  steady  vertical  force  coefficient  if 
there  were  no  reflected  wave  present, fiO   is  the  steady  vertical 
force  coefficient  with  a  reflected  wave  present,  p  is  the 
standing  wave  ratio  given  by 

p  =  HX  -  HM  (A8) 

HX  +  HM 

and  HX  and  HM  are  the  maximum  and  minimum  wave  amplitudes 

that  exist  as  a  result  of  the  summation  of  the  direct  and 

reflected  waves.  HO  used  in  (A2)  is  the  average  wave  amplitude 

given  by 

HO  =  HX  +  HM,  (A9) 

2 

The  complex  potential  functions  given  by  (Al)  and  (A4)  do 

not  satisfy  the  free  surface  conditions  (4)  and  (5) »  hence 

(A3)  and  (A6)  can  be  expected  to  give  reasonable  results  only 

if  ka  «  1  and/or  2kh  »  1,  However,  it  is  assumed  that  (A7)  is 

valid  throughout  the  interest  range  of  ka  and  kh  values „ 
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APPENDIX  B 


Computer  Program  for  Steady  Vertical  Force  Coefficient  Calculation 
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